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Abstract 

We prove that any planar graph on n vertices has less than 0(5.2852") spanning trees. Under 
the restriction that the planar graph is 3-connected and contains no triangle and no quadrilateral 
the number of its spanning trees is less than 0(2.7156"). As a consequence of the latter the grid 
size needed to realize a 3d polytope with integer coordinates can be bounded by 0(147.7"). Our 
observations imply improved upper bounds for related quantities: the number of cycle-free graphs 
in a planar graph is bounded by 0(6.4884"), the number of plane spanning trees on a set of n 
points in the plane is bounded by 0(158.6"), and the number of plane cycle-free graphs on a set 
of n points in the plane is bounded by 0(194.7"). 

1 Introduction 

The number of spanning trees of a connected graph, also considered as the complexity of the graph, 
is an important graph invariant. Its importance largely stems from Kirchhoff's seminal matrix tree 
theorem: The number of spanning trees equals the absolute value of any cofactor of the Laplacian 
matrix of the graph. Furthermore, this number is the order of the Jacobian group of the graph, 
also known as critical group, or as sandpilc model in theoretical physics [21 [3j. This group can be 
represented as a chip firing game on the graph; in this context the number of spanning trees counts the 
number of the stable and recurrent configurations ,4} . The number of spanning trees is also a measure 
for the global reliability of a network. Upper bounds on the number of combinatorial structures 
are often helpful to determine the (exponential) running time for exact algorithms of NP-complete 
problems. 

Our motivation to study the number of spanning trees of planar graphs comes from an application of 
Kirchhoff's matrix tree theorem. Instead of computing the number of spanning trees with Kirchhoff's 
theorem one can use bounds on the number of spanning trees to obtain bounds for the cofactors of 
the Laplacian matrix. These cofactors appear in various settings. For example, Tutte's famous spring 
embedding is computed by solving a linear system that is based on the Laplacian matrix |20|, I21j . As a 
consequence of Cramer's rule the cofactor of the Laplacian matrix is the denominator of all coordinates 
in the embedding. Therefore, by multiplying with the number of spanning trees, we can scale to an 
integer embedding. This idea finds applications in the grid embedding of 3d polytopes [TU [TS] . Before 
we describe this application in more detail, we introduce some notation. 

Let Qn be the set of all planar graphs with n vertices. For a graph G G Qn we denote the number 
of its (labeled) spanning trees with t{G). For every Qn let T(n) be the maximal number of spanning 
trees a graph in this class can have, that is T{n) — max^gg^ {t(G')}. We study the growth rate of the 
function T{n). Since it seems intractable to obtain an exact formula for T(n), we aim at finding a 
value a such that T{n) < a" for n large enough. Notice that the graph that realizes the maximum 
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T{n) has to be a triangulation. Hence, it suffices to look at the subclass of all planar triangulations 
with n vertices instead of considering all graphs in t/„. 

Furthermore, we are interested in the maximal number of spanning trees for planar graphs with 
special facial structure. In particular, we want to bound 

T4(n) — max{t(G)\G is 3-connected and contains no triangle}, 

T^{n) — max {t(G') I G is 3-connected and contains no triangle or quadrilateral}. 

Notice that if a graph is planar and 3-connected its facial structure is uniquely determined |22) . 
Let a" be an upper bound on T4(n) and be an upper bound on T^(n). We refer to the problem of 
bounding a as the general ■problem, and to the problems of bounding 014 and as as restricted problems. 

For embedding 3d polytopes the necessary grid size (ignoring polynomial factors) can be expressed 
in terms of a, 04 and 05. In this scenario we are dealing with 3-connected planar graphs since G is 
the graph of a 3d polytope [18]. If the graph G contains a triangle the grid size is in 0(q;^"), if G 
contains a quadrilateral the grid size is in 0(a4^"). Due to Euler's formula every (3-connected) planar 
graph contains a pentagon - in this case the grid size in 0(a5^"). As a consequence better bounds on 
a, a4 and as directly imply a better bound on the grid size needed to realize a polytope with integer 
coordinates. 

Richter-Gebert used a bound on T{n) to bound the size of the grid embedding of a 3d polytope [15]. 
By applying Hadamard's inequality he showed that the cofactors of the Laplacian matrix of a planar 
graph are less than 6.5". This bound can by easily improved to 6" by noticing that the Laplacian 
matrix is positive semi-definite, which allows the application of the stronger version of Hadamard's 
inequality [IHl page 477]. Both bounds do not rely on the planarity of G, but on the fact that the 
sum of the vertex degrees of G is below 6n. Ribo and Rote improved Richter-Gebert's analysis and 
showed that 5.0295" < T{n) < 5.3" [131 [TB]. The lower bound is realized on a wrapped up triangular 
grid and was obtained by the transfer-matrix method. For the upper bound they count the number 
of the spanning trees on the dual graph. This number coincides with the number of spanning trees 
in the original planar graph. Since the number of spanning trees is maximized by a triangulation, 
the dual graph is 3- regular. Applying a result of McKay |12| . which bounds the number of spanning 
trees on A:-regular graphs, yields the bound of 5.3". Interestingly, this bound is not directly related to 
the planarity of the graph. Therefore, Ribo and Rote tried to improve the bound using the outgoing 
edge approach. The approach involves choosing a partial orientation of the graph and estimating the 
probability of a cycle. To handle dependencies between cycles, Ribo and Rote tried (1) selecting 
an independent subset of cycles, and (2) using Suen's inequality [19]. However, they could only 
prove an upper bound of 5.5202" for T{n), and they showed that their approach is not suitable to 
break the bound of 5.3". For the restricted problems they obtained the bounds T4{n) < 3.529" and 
Ts(n) < 2.847". 

Bounds for the number of spanning trees of general graphs are often expressed in terms of the 
vertex degree sequence of the graph. However, the main difficulty in obtaining good values for a lies 
in the fact that we do not know the degree sequence of the graph in advance. Therefore, these bounds 
are not directly applicable. If we would assume that almost every vertex degree is 6, which is true 
for the best known lower bound example presented in [13] . the bound of Grone and Merris [5] gives 
an upper bound for T{n) of {n/{n — l))"^^6"~^/n, whose asymptotic growth equals the growth rate 
obtained by Hadamard's inequality. To apply the more involved bound of Lyons |llj one has to know 
the probabilities that a simple random walk returns to its start vertex after k steps (for every start 
vertex). Even under the assumption that every vertex has degree 6, it is difficult to express the return 
probabilities in terms of k to obtain an improvement over 6". 

Spanning trees are not the only interesting substructures that can be counted in planar graphs. 
Aichholzer et al. [I] list the known upper bounds for other subgraphs contained in a triangulation: 
Hamiltonian cycles, cycles, perfect matchings, connected graphs and so on. The bounds for Hamilto- 
nian cycles and cycles have been recently improved [5j. 
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Overview. In Section [2] we bound the number of spanning trees by the number of outdegree-one 
graphs, i.e., the number of directed graphs obtained by picking for each vertex one outgoing edge. 
Cycle-free outdegree-one graphs correspond to spanning trees. Therefore we next bound the proba- 
bihty that a random outdegree-one graph has a cycle. For this we analyze the dependencies between 
cycles. In contrast to Ribo and Rote who showed how to avoid the dependencies in the analysis, we 
instead make use of the dependencies. Since our method might also find application in analyzing 
similar dependency structures, we phrase our probabilistic lemma in a more general setting in Sec- 
tion 2.1 More specifically, we develop a framework to analyze a series of events for which dependent 
events are mutually exclusive. In Section [2.2| we apply this framework to bound the probability of the 
occurrence of a cycle. From this we derive in Sections |2.3| and |2.4| a linear program whose objective 
function bounds (the logarithm of) the number of spanning trees. This linear program has infinitely 
many variables, and we instead consider the dual program with infinitely many constraints and present 
a solution in Section [S) 



Results. We improve the upper bounds for the number of spanning trees of planar graphs by showing: 
a < 5.2852, 04 < 3.4162, and as < 2.7156. As a consequence the grid size needed to realize a 3d 
polytope with integer coordinates can now be bounded by 0(147.7") instead of 0(188"). For grid 
embeddings of simplicial 3d polytopes our results yield a small improvement to 0(27.94") over the 
old bound of 0(28.4"). Our improved bound for a implies also an improved bound for the number of 
edge- unfolding cut-trees of a polyhedron of n vertices 7J . 

The maximal number of cycle-free graphs in a triangulation is another interesting quantity. Aich- 
holzer et al. Ij obtained an upper bound of 6.75" for this number. We show in this paper that the 
improved bound for T[n) yields an improved upper bound of 0(6.4948"). 

Multiplying a with the number of maximal number of triangulations a point set can have, gives 
an upper bound for the number of plane spanning trees on a point set. Using 30" as an upper bound 
for the number of triangulations of a point set (obtained by Sharir and Sheffer [17] ) yields an upper 
bound of 0(158.6") for the number of plane spanning trees on a point set. By the same construction 
the number of plane cycle-free graphs can be improved to 0(194.7"). To our knowledge both bounds 
are the currently best known bounds. 



2 Refined outgoing edge approach 

Our results are obtained by the outgoing edge approach and its refinements. For this we consider 
each edge vw of G as a pair of directed arcs v ^ w and w ^ v. Let vi be a designated vertex of 
G, and let V2,...v„ be the remaining vertices. A directed graph is called outdegree-one, if vi has 
no outgoing edge, and every remaining vertex is incident to exactly one outgoing edge. A spanning 
tree can be oriented as outdegree-one graphs by directing its edges towards vi. This interpretation 
associates every spanning tree with exactly one outdegree-one graph. As a consequence the number 
of outdegree-one graphs contained in G exceeds t{G). 

We can obtain all outdegree-one graphs by selecting for every vertex (except vi) an edge as its 
outgoing edge. Let S be such a selection. We denote with di the degree of the vertex Vi. For every 
vertex Vi we have di choices how to select its outgoing edge. This gives us in total rir=2 different 
outdegree-one graphs in G. Due to Euler's formula the average vertex degree is less than 6, and hence 
we have less than 6" outdegree-one graphs of G by the geometric-arithmetic mean inequality. Thus, 
the outgoing edge approach gives the same bound as the strong Hadamard inequality by a very simple 
argument. 

Outdegree-one graphs without cycles are exactly the (oriented) spanning trees of G. To improve 
the bound of 6" we try to remove all graphs with cycles from our counting scheme. Let us now 
consider a random selection S that picks the outgoing edge for every vertex uniformly at random. 
This implies that also the selected outdegree-one graph will be picked uniformly at random. Let P^c 
be the probability that the random graph selected by S contains no cycle. The exact number of 
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spanning trees for any (not necessary planar) graph G is given by 

t{G) = \ \{di\ Pnc- 



2.1 The dependencies of cycles in a random outdegree-one graph 

Assume that the t cycles contained in G are enumerated in some order. Notice that in an outdegree- 
one graph every cycle has to be directed. We consider the two orientations of a cycle with more than 
two vertices as one cycle. Let Ci be the event that the i-th cycle occurs and let Cf be the event that 
the i-th cycle does not occur in a random outdegree-one graph. For events Gi,Cj we denote that 
they are dependent by Ci o Cj and that they are independent by Ci Cj. We say that cycles are 
dependent (independent) if the corresponding events are dependent (independent). 

Two cycles are independent if and only if they do not share a vertex. In turn, cycles that share a 
vertex are not only dependent but mutually exclusive, i.e., they cannot occur both in an outdegree- 
one graph, since this would result in a vertex with two outgoing edges. This gives us the following 
two properties of the events Ci . We say events Ei, . . . , Ei have mutually exclusive dependencies if 
Ei o Ej implies Pr[i?i n Ej] — 0. We say that events Ei,. . . ,Ei have union-closed independencies 
if Ei Ei^, ■ ■ ■ ,Ei ^ Ei^ implies Ei {Ei^ U . . . U Ei^). It is easy to see that the events Q have 
mutually exclusive dependencies and union-closed independencies. 

Lemma 1. // events Ei, . . . ,Ei have mutually exclusive dependencies and union-closed independencies 
then for 1 < k < I 



Pr[f]E^\'f]En<t[ 

j — k i—1 j—k 



( 



V 



Vt[E,] 



n Pr[i?f] n Pr[i?f 

l<i<k: \ k<i<l: 



Before we prove Lemma [T] we prove the following statement. 

Lemma 2. // events Ei, . . . , Ei have mutually exclusive dependencies and union-closed independencies 
then 

l-l 

Pr[Ef I n ^fl < 1 - P^m/ n Pr[Ef]. 



l<i<l: 

E^ ^E, 



Proof. We first prove 



l-l 



Pt[E^ \f]E!]^l- Pr[£;z]/Pr[ f] Ef \ f| Ef]. 



(1) 



l<i<l: l<i<l: 
Ei^Ei Eii^Ei 



The mutually exclusive dependencies imply 

l-l l-l l-l 



V,[Ei I fl Et] = Pr[ii;; n fl E'r]/ Pr[f| E^] = Vi[Ei n fj E^]/ Pr[f| E^]. 



i=l 



i=l 



1=1 



1=1 



The union-closed independencies imply 

Pr[Ein fl Sa=Pr[£;,]Pr[ fj E^], 



= 1,. 



= 1,...,/-!: 

Ei^Ei 
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which by regrouping terms implies equation ([T]). 
A simple consequence of equation ([T]) is that 

i-i 

Pi-[Ef \f]Ef]<l- Pr[Ei] = Pi-[Ef] (2) 

i=l 



It remains to prove that 



pr[ f] E^ n En<PA n ^a- 

l<i<l: l<i<l: l<i<l: 

E^^Ei Ei^E^ E^*^Ei 



Let i' := max{l < i < I \ Ei ^ Ei}. Now, 



Pr[ f] EU n E'r]^Pr[E^,\ fl E^]Pr[ f] Efl fj £f]. 

l<i<l: l<i<l: l<i<l: l<i<i': l<i<l: 

E^^Ei Eii^Ei i^i' E^^Ei E^^Ei 



Thus, by applying Inequality ([2| an induction yields the claim. 

□ 



Proof of LemmaQ After rewriting PrH \ ClU Ef] = Pt[Ei \ ClLl Ef] Pr[n -=1 E^ \ Cltl Ef] 

we obtain by Lemma [2] 



fc-i 



l<i<fc: k<i<j: 

E.^E^ E.^E^ 



The left hand side of equation ([s]) does not depend on the order of the events Ek, ■ ■ ■ Ei. Considering 
once the original order and once the reversed order of these events yields that Pr[n'=fe E^ \ f]^Ii EfY 
equals 

n . P^E,] ' Pv[E,] 



n Pr[i?f] n pr[i?f]-Li n n 

l<i<k: k<i<j: ■> l<i<k: j<i<l- 

Ej^E^ Ej^E^ Ej^E^ Ej^E^ 



Thus, to prove the lemma it suffices to see that for any positive reals a, &, c, (1 — a/b){l — a/c) < 
(1 — a/Vbc)^. By using a/Vbc < [a/b + a/c)/2, i.e., the geometric mean of two positive numbers is 
at most the arithmetic mean, we get, 

2 2 

(1 - a/b){l - a/c) - a/b- a/c<\+°- - 2a/\/bc = (1 - a/\/bcf. 

be be 

□ 



2.2 Bounding the probability of the appearance of cycles 

Before estimating the probability Pnc in terms of the vertex degrees, we introduce some notation. A 
cycle of length k is called a k- cycle. The k- extension of a cycle is the union of a cycle with all its 
dependent /c-cycles. We say that the degree of a cycle is the ordered sequence of the degrees of its 
vertices. Let Cabc be a 3-cycle spanned by Va,Vb, Vc, and let the degree of Cabc be {da, di,, dc) = («, j, k). 
We denote the degrees of the vertices adjacent to Va that are not part of Cabc by the sequence A. In 
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the same fashion we denote the degrees of the vertices around Vb by B and the around Vc by C. The 
ordering in A, B, C respects the counter clockwise ordering of the vertices around Va^ Vb, Vc in a planar 
embedding. Since G is planar and 3-connected the ordering of the sequences is uniquely determined 
up to a global reflection [22 . Notice that a vertex might occur in two different sequences. We call the 
tuple (i, J, fc, A, B, C), the signature of the 2-extension of Cabc- Similarly, we define the signature of a 
2-extension of a 2-cycle Cab by the tuple A, B). The naming convention is depicted in Figure[T] 

(13 — bi 

ai = On \ / 62 / \ 



/j = 6 



64 = ci , 



..fc = 7 



(5,6, (ai,..., 04), (&!,..., 65)) 



(5, 6, 7, (ai, . . . , as), (61, ... , 64), (ci, . . . , C5)) 



Figure 1: Convention for naming the signatures of 2-extensions of 2-cycles (on the left) and 3-cycles 
(on the right). 



We can express P^c as Pr[nj=i ^'f]- Our goal is to apply Lemma jl] to bound this probability. As 
a first step we discuss how to express the number of spanning trees t{G) when the distribution of 
signatures of G is known. Instead of t{G) we bound its logarithm, i.e., 

n t 

logi(G)=^logd,+logPr[f|Cf]. (4) 

i=2 ]=1 

The probability that an event Ci occurs can be expressed in terms of vertex degrees. In particular, 

Pr[Ci] — l/{dadb) the i-th cycle is a 2-cycle on the vertices tiatife, 

Pr[Cj] = 2/ Jl^. ,y g^c?a the j-th cycle is at least a 3-cycle on the set Z. 

The way we proceed depends on whether we are addressing the general problem (i.e., we want to 
bound a) or one of the restricted problems (i.e., we want to bound 04 or a^). In the latter case we 
limit our analysis to 2-cycles only. In the general case we consider all cycles of length 2 and cycles of 
length 3 that are triangles in G. 

We start with the general problem. Assume that all cycles C are enumerated such that the first 
ts cycles are the triangles in G, and the last ^2 cycles are the 2-cycles of G. In total we consider 
t := t2 + ^3 cycles. All remaining cycles are ignored. Discarding the larger cycles gives an upper 
bound on P^c and is therefore applicable. We apply Lemma [l] with k — 1 and I — to bound 
P^injLi C'j]; which is the probability that no 3-cycle occurs. To take also the 2-cycles into account 
we consider the probability that no 2-cycle occurs under the condition that no triangle occurred as 3- 
cycle, which is Pr[nj=t3+i Cjl Hj-Li C^]- Notice that this probability has the form stated in Lemmajl] 
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for I = t and fc = ^3 + 1. Thus, we can bound log Pr[p|*^-^ C|] from above by 

PS V n Fv[cnj jir+i v n priQi n priQiy 

Wl<i<t3: l<i<t3 + l: \ t3<i<t: 

Equation ([5]) is a sum over cycles. Each summand in this sum depends only on the signature of the 2- 
extension of such cycle. Hence, we can group the summands in ^ with identical signatures. We denote 
the number of 2-extensions of 2-cycles with signature {i,j,A,B) by the variable fij{A,B). Similarly, 
the number of 2-extensions of 3-cycles with signature {i,j,k,A,B,C) is denoted by fijk{A, B,C). 
In order to simplify matters, we refer to fij{A,B) and fijk{A, B,C) simply as fij and fijk, or as / 
variables. 

For better readability we introduce the following notations {X is used as a placeholder for A, B, 
or C, and X clS 3t placeholder for a, b, or c): 

P2ir,X) n fl-— )' ^3(r,X):= [] — 

l<p<r-l ^ ^^P^ l<p<r-2 ^ T-XpXp+i 
Pjj(A,B) ;= 1 ^ r . 

^jp^ii^^ A)p,{h B)(i - ^) (1 - j^j ,/m:A)P2{j;B) 

P,,k{A,B,C) := 1- 



ijk^P^ii, A)P,ij, B)P^ik, 0(1- ^)( 



1_ 2 ^_ 2 



We rephrase ([5]) as 

t 

logPr[f| q] < E B, C) logPyfc(A, 5, C) f,,{A, B) \ogP,,{A, B). (6) 

j = l i,jM,A,B,C i:j,A,B 

The sums in the last expression (and following similar sums) range over all feasible signatures. Let us 
now consider the restricted problems. Both restricted problems are easier to analyze than the general 
problem, since we consider only 2-cycles. To bound PriPlj-Li &Pply Lemma jlj with k — 1 and 

I = t2- Following the presentation of the general problem we define 

P^j{A,B) 1- 



i^P2{i,A)P2{j,By 



and obtain for the restricted problems 

logPr[f]q]< hiA,B)logh,{A,B). (7) 

j=l id,A,B 



2.3 A charging scheme for the vertex degrees 

If we insert the bounds ([6| or Q into equation Q we obtain an upper bound for t{G) in terms 
of the signatures of G. However, we would like to express the first part of equation Q, which is 
D := X]"=ilogdi, also in terms of the / variables. For convenience we include logdi in the sum for 
D, which is applicable since we are looking for an upper bound. 

Let us first discuss the general problem. We split D into four parts: Di := HiD for i = 1, . . . , 4, 
with Mi — 1- The parameters /i^ will be fixed later. We express Di and D2 by the fij variables 
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and D3 and D4 by the fijk variables. Every vertex Va contributes fii log do to Di. On the other hand, 
every vertex Va is part of da 2-cycles. We charge the total amount of /ii log da uniformly to these 
2-cycles. Thus, every 2-cycle incident to Va gets fiilogda/da from Va- In a similar fashion we charge 
D2 to the 2-extension of 2-cycles. Let VaVi, be an edge in G and let Vj. 7^ Wf, be a vertex adjacent 
to Va- Distributing fi2^ogdr uniformly, assigns every 2-extension with "endpoint" Vr the fraction of 
fi2\ogdr/{dr{da — 1)) from Vr- For and D4 we argue analogously. We can therefore express D by 

4ogi log.r 



J 



i,j,A,B 

\ - ^ ^ ^ f \ - log ttr \ - log br 



. 0,r{i — 1) ■'^^ brii — 1) 



logz ^ logj ^ logfc 

i 7 fc 

i,j,k,A,B,C ■' 



(8) 



We can now express log Pnc as sum over all signatures. This sum can be subdivided into one part 
that contains the fij variables and one part that contains the fijk variables. The part that considers 
the 2-cycles is given by 

D1 + D2+ J2 .nj{A,B)logP,,{A,B), (G2) 

i,j,A,B 

and the part that considers the 3-cycles is given by 

D5 + D4+ njk{A,B,C)logP,jkiAB,C). (G3) 

i,j,k,A.B,C 



For the restricted problems we only have 2-cycles. Using bound ([t]) and setting — iJ-i = 0, we can 
bound the number of spanning trees by 

D1 + D2+ J2 h{A,B)logP,,(A,B). (R2) 

i,j.A,B 

2.4 Finding constraints 

In this section we construct necessary conditions for the / variables that have to hold for planar 
graphs with n vertices. We reuse the ideas from the charging scheme in Section [2?3l Instead of giving 
every vertex logdj to distribute, we assign to every vertex an amount of 1. This gives us a total of n 
units. Following the construction of the equations of ([s]) we obtain 

Y h{A,B)[-^+^) =n, (A2) 

i,j,A,B ■' 

Y hk{A,B,C){^ + - + l) ^n, (A3) 

i,j,k,A,B,C * 
i,j,A,B a^eA '^^ ' b^eB ' 
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Another set of constraints is given by the number of 2-cycles and 3-cycles a planar graph can have, 
which is related to the number of edges and faces of G. Every 2-cycle is counted by some fij variable, 
hence the sum over all fij equals the number of edges, which we name m. Since we consider only 
3-cycles of triangles, the sum of the fijk variables equals the number of triangles, which for a planar 
graph is at most 2n. We obtain 

E hiAB)<m, (C2) 
hk{A,B,C)<2n. (C3) 

i,j,k,A,B,C 



For the general case we have m < 3n, for the restricted case where quadrilaterals are allowed we have 
m < 2n, and for the remaining case we have m < 5n/3. All these bounds can be obtained by a simple 
double counting argument using Euler's formula. As trivial condition we restrict the / variables to 
be non-negative. 

The constraints so far might be fulfilled by a signature that does not come from a planar graph. 
In particular, the degree sequence of the graph induced by the cycles might be unrelated to the degree 
sequence of the graph induced by the 2-extensions. To overcome this ambiguity we consider the 
number of edges with vertex degree i at one vertex and degree j at the other. Let this number be 
Uij. Clearly, we have that riij = ^ fij {A, B), where the sum ranges about all feasible sequences 
A, B. On the other hand, can be counted by its appearances in the 2-extensions of 2-cycles. Every 
edge with degree (?, j) will show up in (i — 1) -I- {j — 1) 2-extensions. Let Xii^) denote the number of 
appearances of i in the sequence X. We can express ((z — 1) -f (j — as J2k A b fik{A, B)xj{A) + 

A B fkj{A, B)xi{B). This leads us to a new constraint of the form 

(z-f j-2)^/,,(AS)= E lMB)xj{A)+ fkAA,B)xr{B). (Ezj) 

A,B k,A,B k,A,B 

In the case where the smallest face of the graph is a pentagon, we were able to improve the solution 



of the linear program by adding the constraint (E33) . Other constraints of the form ( Eij ) gave no 
improvement. 



3 Solving the linear programs 

For the general case and a fixed n we can solve the problem of maximizing t{G) while fulfilling the 
constraints of Section |2.4| by solving two linear programs. The first problem, which is derived from 



the 2-cycles, is given by the objective function (G2|, and the constraints (A2-C2). The problem 



associated with the 3-cycles consists of the objective function (G3|, and the constraints (A3-C3). Our 
goal, however, is not to solve these problems for a fixed n, but to study the solution in terms of n. 
To get an expression in terms of n, we normalize the / variables. Instead of letting / count absolute 
numbers, we consider / as ratio between these absolute numbers and n. This allows us to cancel n in 
every constraint and objective function. 

By removing n we lost a natural bound on the largest vertex degree, which is rt — 1. Moreover, we 
cannot bound the number of signatures and have to deal with infinitely many variables. To overcome 
this difficulty we look at the dual program. The dual has up to four variables, which we name Ai,A2,A3, 
and A4. Since we cannot bound the number of variables in the primal, we have infinitely many dual 
constraints. By the weak linear programming duality it suffices to find a point in the feasible area of 
the dual, because every such point is an upper bound for the primal solution. Appendix A lists all 
relevant dual programs. Let Z2 be the LP solution for (G2) and Z3 be the LP solution for (G3). The 
combined bound for a can be computed by adding exp(Z2) and exp^Z^). 
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LP instance 


/I values 


solution 


a 


genercil problem (2-cycles) 
(G2)^ max, s.t. (A2-C2), m = 3 


Ml = 0.3, M2 = 0.25 


A3 = 0.232445 


5.28515 


general problem (3-cycles) 
fG3W max s t fA3-C3) 


Ma = 0.225, M4 = 0.225 


Ai = 0.0980332 
A2 = 0.192612 
A3 = 0.247984 


restricted problem (no triangles) 
(R2)^ max, s.t. (A2-C2), m = 2 


Ml =0.9,M2 =0.1 


A3 = 0.614264 


3.41619 


restricted problem (no 3,4-gons) 
(R2)^ max, s.t. (A2-C2,E33) 
m = 5/3 


Ml = 0.615, M2 = 0.385 


Ai = -0.615054 
A3 = 0.744706 
A4 = 0.001954 


2.71567 



Table 1: The results of the dual linear programs. Not specified A values are zero. 



We use a two step approach to find the solution of the dual programs. First we solve the dual 
program with a finite number of the constraints. We expect that the maximum number of spanning 
trees will be realized on a graph with evenly distributed vertex degrees. Therefore, we expect in the 
general case that the variables fij, with i, j far away from 6 will be zero, and hence the corresponding 
dual constraints will not hold with equality. For the restricted problems we expect values for i,j 
around 3 and 4. We use these assumptions to construct (finite) linear programs that will most likely 
contain the constraints that determine the dual solution. We prove later that the other constraints 
are also fulfilled. The candidates for the dual solution that were obtained by the finite programs are 
listed in Table |3] The table also lists our choice of m values. We picked these particular numbers 
because they have been experimentally proven to be useful in the later analysis. 

The verifications of the LP solutions is tedious. We redirect the interested reader for a com- 
plete analysis to the Mathematica scripts, downloadable under http://wwwmath.uni-muenster.de/ 



u/schulaxid/scripts-nst .zip For the 2-cycle part of the general problem we present the general 
idea behind the analysis in this place. Remarks for the analysis of the other programs can be found 
in Appendix B. 

Lemma 3. For mi = 0.3, M2 = 0.25 and Ai = 0, A2 = 0.180948, A3 = 0.232445, and m ^ 3 all dual 
constraints (A2-C2) for feasible signatures are fulfilled. 

Proof outline. The complete proof including technical details and algorithms for brute force testing 
of some cases can be found in the Mathematica scripts. 



We have to show that the left-hand side of ( 10 ) is negative. Let us ignore the log-term for now 



it is negative, and therefore it suffices to show that ( 10 ) without the log-term is negative. Let a 



be some vertex degree stored in the sequences A or B. We notice that the "effect" of a in (10) is 



M2 log(a)/(a(i — 1)) — A2/(a(i — 1)). For integers this term is maximized for a = 6, and for a < 5 



monotonically increasing. Hence, it is sufficient to show that (10 1 holds for all a > 6. 
Using the lower bound a = 6 gives us 



VeA ^ ^ beB 'J \a<^A ^ ' beB ^-^ ' , 



< 2M2(log6)/6-2A2/6. (9) 



We know that Ai = 0, and that (logj)/j is maximized for integers at j — 3. Hence, if the expression 
Mi((logi)/i -f (log3)/3) -I- 2M2(log6)/6 — 2A2/6 — A3 is negative, the dual constraint is fulfilled. By 
considering the partial derivative in j, we see that this holds for all z > 31 (and a symmetric argument 
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implies the same for j > 31). For all other pairs ij we check the dual constraint by hand using again 
the lower bound 6 for the entries of A and B. By this we can eliminate more cases and are left with 
83 tuples (see Appendix B). 



We will now consider the log-term of ( 10 1 again. Therefore, we have to obtain new upper and 
lower bounds for the entries of A and B. Let a be a vertex degree that is stored in A or i?. It can 
be observed that the log-term in terms of a is increasing. Since the function without log-term was 
increasing between 3 and 6 the function with log-term is increasing as well in this range and hence 
the maximum has to be attained on a value larger than 6. 

To obtain an upper bound on the entries of A and B we argue as follows: We first assume that all 



entries of A and B equal 6, except for some & A. We look at the derivative of the constraint (10) 
in Qx and observe that the maximum is realized at some value a^. We study now what happens if we 
change one of the values, say ay, which we assumed to be 6. Let g{ax,ay) be the dual constraint in 
ay and a^- It can be shown that g{ax, ay) / duxday is negative. Hence, for ay > 6, 

Q ^ dg{ax,6) _ dg{ax,ay) ^ {dg(ax,6) - g{ax,ay)) 
dax dax 

In other words, the differences between g(ax,6) and g(ax,ay) are increasing in ax- Let a'y be 
fixed and d = g(a^, 6) — g{a'^,ay). Clearly the function h{ax,a'y) := gia^, a'y) + d attains its extrema 
on the same positions as g{ax^a'y) and its differences to g{ax,Q) are also increasing in ax- Since 
h{a'^,a'y) = g(a^,6), the graph of h{ax-,a'y) is always below g{ax,Q) for ax > a'^- Therefore, no 
maximum can be attained on h{ax,a'y) for ax > a'^ and hence the smallest upper bound is achieved 
for ay = 6. Hence a'x is a valid upper bound for ax- The computed bounds for every remaining pair 
ij can be found in Appendix B. 

As last step we check all possible signatures that are left to test. Notice that there are still many 
cases open. We observe that the ordering of the degrees in A and B does not matter, except for 
P^{i,A)P'j,{j,B)(^ — jj^^ ^1 — jjbl^- But this expression can be bounded by using the maximum 
in A for all ax, and the maximum in B for all bx- With this estimation we can try all sequences 
with sorted sequences A and B- This eliminates most of the cases. There are three pairs left to test. 



namely (6, 6), (5, 6), and (5, 7). We test the constraints for the remaining signatures by checking ( 10 ) 
for all combinations. 

The solutions for the linear programs lead to the main theorem. 

Theorem 1. Let G he a planar graph with n vertices. The number of spanning trees of G is at most 
0(5.28515"). If G is 3-connected and contains no triangle, then the number of its spanning trees is 
bounded by 0(3.41619"). If G is 3-connected and contains no triangle and no quadrilateral, then the 
number of its spanning trees is bounded by 0(2.71567"). 



4 Further bounds and future work 

The results of Theorem [T] improve several related upper bounds. Using the observations by Ribo et 
al. [13] we obtain the following bounds for grid embeddings of 3d polytopes. 

Corollary 1. Let G be the graph of a 3-polytope V with n vertices. V admits a realization as 
combinatorial equivalent polytope with integer coordinates and 

L no coordinate larger than 0(147.7"), 

2. no coordinate larger than 0(39.9"), if G contains a quadrilateral, 

3. no coordinate larger than 0(28.4"), if G contains a triangle. 
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The number F(n) of cycle-free graphs in a planar graph with n vertices is bounded by the number 
of selections of at most n — 1 edges from the graphs [T]. Thus, F{n) < J2k=o < g < 1/2 

we have that ^[^"J (J^) < 2^(«)", where H{q) := -log(g)« - log(l - is the binary entropy 

(see for example [H page 427]). This shows that, F{n) < 6.75" by setting m — in and q — 1/3. 

We give a better bound based on the bound for the number of spanning trees. We first bound 
the number k) of forests in f/„ with k edges. On one hand, the above argument yields an 
upper bound of F{n,k) < /i(fc) := (^\7®). On the other hand, every forest with k edges can 
be constructed by selecting k edges from a spanning tree of C/„. This gives as alternative bound 
F(n, k) < f2{k) —: ("^^)T(n). Now, the number of cycle-free graphs is bounded by 

n-l 

F{n) = y F{n,k) < n max F{n,k) < n max min(/i(gn), /2(gn)). 

— ^ 0<fc<n 0<o<l 

fe=0 

We use (^'^) < 2"^^'^^ as upper bound for the binomial coefficients (see for example fSJ page 1097]) to 
obtain 

/i(9ri)</i(gn):=23"«('/3) and h{qn) < h{qn) -.^ T{n)2^"('^l 

The computed maximal value for the minimum of /i and /2 is realized at qn = 0.94741 n. This yields 
a bound of n ■ 6.4948" for the number of cycle-free graphs. For the computation of these values we 
used numerical methods. Observe that fi{qn) realizes 6.4948" at a larger value q than f2{qn). The 
correctness of the numerical computations follows from the monotonicity of fi and /2 in (n/2, n]. For 
the number of plane spanning trees and cycle- free graphs on a planar point set, we obtain improved 
upper bounds by multiplying our bounds with the bound of 0(30") on the maximum number of 
triangulations on a planar point set |17j . 

Theorem 2. The number of cycle- free graphs in a planar graph with n vertices is bounded by n ■ 
6.4948". The number of plane spanning trees on n planar points is in 0(158.6"), the number of 
cycle-free graphs in 0(194.7"). 

We expect better bounds for the number of cycle-free graphs in a planar graph from a more direct 
application of the outgoing edge approach. By adding a new vertex that is linked to a subset of the 
other vertices, every cycle-free graph can be turned into a spanning tree of the augmented graph. 
Without excluding any cycles we get a bound of 7". Under the assumption that almost every vertex 
has degree 6, the refined outgoing edge method would yield a bound of 6.5027" when excluding 2- 
cycles and of 6.4244 when excluding 2 and 3-cycles. So far we were not able to check all constraints 
of the corresponding linear programs. 

We finish our presentation with a discussion on how one could improve our results further. Since 
we consider only 2-cycles and 3-cycles from triangles, one would obtain a better bound for P^c by 
taking also larger cycles into account. We do not expect to win anything by considering 3-cycles that 
are not triangles, because in the lower bound example (the wrapped up triangular grid) all 3-cycles are 
triangles. The analysis using larger cycles is more complicated and needs an extensive case distinction. 
Furthermore, we expect that there would be too many cases left for the brute force check. From our 
perspective, the following refinement seems tractable: Beside the 2-cycles, and 3-cycles on triangles, 
we also analyze 4-cycles that belong to two triangles sharing an edge (the diagonal) . The 4-cycles can 
be analyzed by extending the events Ci for the 2-cycles to the following event: the i-th 2-cycle occurs, 
or the corresponding 4-cycle, whose diagonal is associated with the i-th cycle occurs. Assuming that 
the solution of the corresponding linear program is given by having almost every vertex degree 6, this 
would lead to a = 5.25603. Since the resulting linear program is more complicated, the verification 
of the dual solutions is tedious. 

Notice that Lemma [T] uses two enumerations of the events Ci to avoid the infiuence of the ordering. 
An elaborated enumeration scheme of the events Ci might give better bounds. Furthermore, we could 
consider "extension of extensions" to analyze larger locally connected pieces of the graph at once. 
This results in a powerful but very complicated incarnation of the outgoing edge approach. 
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The reader might think, that additional constraints in the hnear programs might improve the 
outcome of our analysis. However, we expect that the solutions of the dual programs give the correct 
distribution of signatures. In particular, the solutions the dual programs match the candidates for the 
lower bound examples that were presented in jl3| . 

Acknowledgements: We thank Giinter Rote for suggesting this problem to us and for many inspiring 
and fruitful discussions on this subject. 
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Appendix A: Dual programs 

General Problem (2-cycle part): 

Minimize Ai + A2 + 3A3, 
such that, A3 > 0, and for all signatures {i,j,A,B): 



log {A, i?) + Ml (— + — j + ( ^ + :L ^y^I) j 

General Problem (3-cycle part): 



Minimize Ai + A2 + 2A3, 
such that, A3 > 0, and for all signatures (i, j, k, A, B,C): 



Restricted Problem (no triangles): 

Minimize Ai + A2 + 2A3, 
such that, A3 > 0, and for all signatures {i,j,A,B): 

'logz logj\ / >r-^ log a,. ■s-^ log6r 



Restricted Problem (no triangles, no quadrilaterals): 

Minimize Ai + A2 + 5/3 A3, 
such that, A3 > 0, and for all signatures A, B): 



1 i> lA T}\^ /logi logj\ / loga^ logfer \ 

iogP,,(AS) + m(— + — j+M2(l^ ^;^+ L lyj 



where Eij{A, B) — < 



if « > 4, 

-X3(^) ifz = 3,j>4, 

[4-X3(A)-X3(P) ifj = 3,j = 3. 



15 



Appendix B: Verification of the solutions of the dual linear 



programs 

We give in this place informations how the candidates for the dual solutions in Table [3] can be 
checked. The complete analysis was computed by a series of Mathematica scripts, which can be found 
http://wwwmath.uni-muenster.de/u/schuland/scripts-nst.zip. These scripts also perform the 
necessary brute force part of our analysis. We follow the general, more detailed, procedure of verifying 
the dual solutions presented in the proof of Lemma |3] We discuss the four linear programs one by 
one. 

General case (2-cycle part): This case is covered in detail in the proof of Lemma [sj We only list in 
Table [2] the intermediate results of the analysis, i.e., the cases which need to be checked brute force. 
The table contains also the upper bound for the entries of A and B for each case. Notice that the 
calculations for the bounds of the first entry oi of A and of the first entry bi of B are slightly different. 
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Table 2: Cases for the brute force test in the 2-cycle part of the general problem. 
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General case (Z- cycle part): The lower bound for the entries of the entries oi A,B,C can be obtained 
by maximizing (/i2loga — A2)/a, which is maximized for integers at a = 6. To get an upper bound 
for the j, k values compute first the maximum of /Lti \ogx/x — Xi/x —: S. Knowing S and the lower 
bound for for the entries of A, B, C, we are able to compute the the upper bound for the k values 
as the maximum of ^i(logi/i + 2 log S/S) + 3/i2 log 6/6 — Xi{l/i + 2/d) — 3A2/6 — A3. This expression 
is negative for all integers gre ater than 13, hence we have to test all remaining tuples for k < 13. 



These cases arc listed in (11) as set of triplets i,j,k. For all i,j,k triplets we obtained that 7 is an 



upper bound for the entries in A, B, C. 

{(4, 4, 4), (4, 4, 5), (4, 5, 5), (5, 5, 5), (4, 4, 6), (4, 5, 6), (5, 5, 6), (4, 6, 6), (5, 6, 6), (6, 6, 6), 

(4, 4, 7), (4, 5, 7), (5, 5, 7), (4, 6, 7), (5, 6, 7), (6, 6, 7), (4, 7, 7), (5, 7, 7), (6, 7, 7), (7, 7, 7), 
(4, 4, 8), (4, 5, 8), (5, 5, 8), (4, 6, 8), (5, 6, 8), (6, 6, 8), (4, 7, 8), (5, 7, 8), (6, 7, 8), (7, 7, 8), 
(4, 8, 8), (5, 8, 8), (6, 8, 8), (4, 4, 9), (4, 5, 9), (5, 5, 9), (4, 6, 9), (5, 6, 9), (6, 6, 9), (4, 7, 9), 

(5, 7, 9), (6, 7, 9), (4, 8, 9), (5, 8, 9), (4, 4, 10), (4, 5, 10), (5, 5, 10), (4, 6, 10), (5, 6, 10), 
(6, 6, 10), (4, 7, 10), (5, 7, 10), (4, 4, 11), (4, 5, 11), (5, 5, 11), (4, 6, 11), (5, 6, 11), (4, 4, 12), 

(4, 5, 12), (5, 5, 12), (4, 6, 12), (4, 4, 13), (4, 5, 13), (3, 4, 5), (3, 5, 5), (3, 4, 6), (3, 5, 6), 
(3, 6, 6), (3, 4, 7), (3, 5, 7), (3, 6, 7), (3, 7, 7), (3, 4, 8), (3, 5, 8), (3, 6, 8), (3, 7, 8), (3, 4, 9), 

(3, 5, 9), (3, 6, 9), (3, 4, 10), (3, 5, 10), (3, 4, 11)} (11) 



Restricted cases: For the restricted cases we use the natural lower bound of 3 for the entries of A and 
B. If the smallest face of G is a quadrilateral, an upper bound for i,j can be computed by maximizing 
^i((logi)/z + log(3)/3) + 2^2(log3)/3 — A3. This expression is maximized for integers when i — 10, and 
hence negative for i > 10. If the smallest face of G is a pentagon, we have to compute the maximum 
of fii{\ogx)/x — Xi/x =: S first. Using S for j we can show that ^i(logi)/i + (log5)/J) + 2^2(log3)/3 — 
Al(l/^ + 1/6) — 2A2/3 — A3 is negative for all j > 11 . Thus were are left with trying out all cases for 
i,j < 10 (if G contains no triangle) and i,j < 11 (if G contains no triangle and quadrilateral). The 
upper bounds of the entries of A and B are computed for each of these cases as explained in the proof 
of Lemma [3l 

The "pentagon case" contains one more subtlety. The constraints where i or j are 3 are slightly 
more complicated because we included (E33) in the primal program. For convenience we assume 
that the upper bound on the entries of A and B for these ij values is at least 4. This simplifies 
the computations since for a distinct neighbor the influence of (E33) is a constant when varies 
between 4 and 00. See Table |3] for the remaining cases. Notice that we can assume that the sequences 
A and B are ordered since the order of A and B is not relevant for the dual constraints of the restricted 
cases. 
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test of the restricted problems. 
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